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Abstract: We carry out a mathematically rigorous investigation into the equi- 
librium thermodynamics of massless and massive bosons confined in general- 
ized Sicrpinski carpets (GSCs), a class of infinitely ramified fractals having non- 
integer Hausdorff dimensions dh- Due to the anomalous walk dimension dyj > 2 
associated with Brownian motion on GSCs, all extensive thermodynamic quan- 
tities are shown to scale with the spectral volume with dimension dg = 2{dyi/d^) 
rather than the Hausdorff volume. We prove that for a low-temperature, high- 
density ideal massive Bose gas in an unbounded GSC, Bose-Einstein condensa- 
tion occurs if and only if dg > 2, or equivalently, if the Brownian motion on the 
GSC is transient. We also derive explicit expressions for the energy of blackbody 
radiation in a GSC, as well as the Casimir pressure on the parallel plate of a 
fractal waveguide modelled after a GSC. Our proofs involve extensive use of the 
spectral zeta function, obtained via a sharp estimate of the heat kernel trace. 
We believe that our results can be verified through photonic and cold atomic 
experiments on fractal structures. 
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1. Introduction 

This paper is devoted to the study of the thermodynamics of quantum gases in 
fractal spaces. It was long recognized, by H. Lorentz and H. Weyl |31j, that a 
deep connection exists between the thermodynamic properties of quantum gases 
and the underlying spectral geometry. One can probe the asymptotic behavior 
of elliptic linear differential operators in a given space by measuring the energy 
or pressure of a quantum gas in the said space. Since then much of this work has 
been done on Euclidean spaces. Here we wish to address the following question: 

How does the fractal geometry affect the laws of quantum many-body physics? 
Conversely, what information about the fractal geometry can we obtain using 
probes made of quantum particles? 

Below are some specific examples which are familiar to students of quantum 
and statistical physics, and which have wide ramifications in condensed matter 
physics and quantum field theory. 

— Does Bose-Einstein condensation of atomic gas occur in non-integer dimen- 
sion? If so, what dimension is it? And what is the critical density for conden- 
sation? 

— What is the analog of the Stefan-Boltzmann law for blackbody radiation when 
the blackbody itself is a fractal? 

— At zero temperature, what is the Casimir pressure produced by vacuum fluc- 
tuations inside a fractal? 

We will answer all three questions in the case where the fractal is a generalized 
Sierpinski carpet ( GSC), whose representatives are the standard two-dimensional 
Sierpinski carpet and the three-dimensional Menger sponge (Fig. [l]) . These frac- 
tals have connected interior, are highly symmetric, and most of alTfare infinitely 
ramified, which makes the analysis difficult. Conventional analytic tools, such 
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as Fourier transform or spectral decimation, no longer apply. Essentially all the 
rigorous results known today originate from the study of Brownian motion. On 
the other hand, GSCs are embeddable in Euclidean space and resemble more 
realistic fractals found in nature. So we believe that a careful analysis of quan- 
tum gases in GSCs is warranted, because it provides us an avenue of attacking 
quantum many-body problems on general irregular spaces. 

Our exposition is aimed at both physicists and mathematicians. It serves 
a dual purpose: to highlight the latest developments from the mathematical 
analysis on fractals, which have just begun to percolate through the physics 
community; and to illustrate how state-of-the-art potential theoretic results can 
be applied to answer physically inspired problems. 

Notations. We write C and c for positive constants which may change from 
line to line. If a constant has a specific value, then we will add a numeral subscript 
to indicate this, e.g. Ci. Given two real- valued functions / and g, we say that 
/ is comparable to g if there exist constants c, C > such that cg{x) < f{x) < 
Cg{x), denoted by f{x) x g{x) for short. For any A CM.'^ and L > 0, we write 
LA = {x eR'' : x/L (E A}. Finally, we use Ng and N to denote, respectively, the 
set of natural numbers with and without 0. 



1.1. Generalized Sierpinski carpet. Let Fq := [0,1]"^ be the unit cube in M"^, 
d > 2. Fix a length scale factor Ip G N, Zi? > 3, and let iS„ be the collection of 
closed cubes of side Ip"' with vertices in lp"Z'^. For A c M'', let 5„(A) = {5" e 
Sn '. S G A}. Denote by the orientation-preserving affine map which maps 
Fo to 5 G Sn. 

Introduce a decreasing sequence {F„}„ of closed subsets of Fq as follows. Fix 
mp G N, 1 < mp < Ip, and let Fi be the union of m^^ distinct elements of 
Si{Fq). In other words, we construct Fi by retaining m^? of the cells of length 
Ip , and removing the rest. We will refer to mp as the mass scale factor. Then 
by iteration we let 



Fn+i= U -fs(i^i)= U 1's{Fn),n>l. 

ses„{F„) seSi{Fi) 
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This iterated function system of contractions {^s} has a unique fixed point 
F = standard arguments, the Hausdorff dimension of F is dii{F) — 

logmp/^oglF- 

Definition 1.1. We say that F is a generahzed Sierpinski carpet (GSC) if the 
following conditions on Fi hold: 
(HI) (Symmetry) Fi is preserved under the isometries of the unit cube. 
(H2) (Connectedness) Int(i^i) is connected, and contains a path connecting the 

hyperplanes {xi = 0} and {xi = 1}. 
(H3) (Non-diagonality) Let B be a cube in Fq which is the union of 2'* distinct 

elements of Si. Then «/Int(Fi n B) 7^ 0, it is connected. 
(H4) (Borders included) Fi contains the segment {x G M.'^ : Xi G [0, 1], X2 = ■ ■ ■ = 
xd = 0}. 

We will denote by dF the boundary of F, and by doF — d{[0, 1]'^) n dF the 
outer boundary. 

In order to discuss the thermodynamic limit we need to consider unbounded 
versions of GSC. For each GSC F, we call F :— IJ^g IpFn the corresponding 
pre-carpet, and F^o ■= U^o unbounded carpet. Plainly speaking, the 

former is the infinite "blow-up" of the fractal, whereas the latter has both "blow- 
up" and "blow-down." Our focus will be on the unbounded carpet, where sharp 
results can be stated, though we believe many of them carry over to the pre- 
carpet modulo minor modifications. 



1.2. Statement of the problem. Let be a GSC, be the self-similar measure 
on F, and Hi := L'^{F^v) be the single-particle Hilbert space. Given a self- 
adjoint, bounded-below Hamiltonian operator H : Hi — > Hi, which captures 
the physics of a certain ideal Bose gas inside F, we wish to compute its grand 
canonical partition function ^p^p, at inverse temperature /3 e (0, 00] and chemical 
potential jj, € (— 00, inf Spec(i/)]: 

log H^,^ = -Trm log (1 - e-^^"-^A . (1.1) 



Note that (1.1) already takes into account the Bose-Einstein statistics satisfied 
by the particles (see Section |4]). 

The Hamiltonian considered in this paper involves the (nonnegative) Lapla- 
cian —A on the measure space {F,v), constructed via either the Barlow-Bass 
approach (Brownian motions on the outer approximations) or the Kusuoka-Zhou 
approach (random walks on the associated graphs). Up to time change, both 
Laplacians generate the same (and the unique) Brownian motion which is in- 
variant under the isometries of the carpet, as |'6! has shown. We will concentrate 
on the following quantum gases, working in units where Planck's constant (K), 
and either the speed of light or twice the atomic mass, are 1. 

— Massive Bose gas {e.g. atoms), where H — —A and /j, < inf Spec(iJ). The 
associated scalar field ■0 satisfies the usual Schrodinger equation idtip = —Aip. 

— Massless Bose gas {e.g. photons), where H = \J —A and /i = 0. Recall 
that the massless scalar field i\) satisfies the wave equation —dtt'4' — {—A)ip. 
Since the Laplacian —A is a nonnegative self-adjoint operator, it has a unique 
square root \/—A, also known as the Dirac operator. Thus we may take the 
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formal square root of the wave equation to obtain a Ist-order-in-time PDE 
idtip — yJ—Aip, which is equivalent to a Schrodinger equation idtil> — Hijj 
with H = 



1.3. Methodology. An immediate obstacle to a rigorous thermodynamic calcula- 
tion on GSCs is the lack of knowledge about the exact spectrum of the Lapla- 
cian. In fact, even the eigenvalue counting function (or the integrated density 
of states in physics parlance) in its optimal form cannot be proved (see Re- 



mark 2.3). Fortunately, many efforts which went into studying the short-time 



asymptotics of the heat kernel on GSC have produced sharp enough estimates, 
to the point that we can show that the corresponding spectral zeta function 
C,a{sj"i) = Tr^^(— Z\ -I- 7)"" admits a meromorphic extension (Section [3|. This 
zeta function is then used effectively in computing the thermodynamic partition 
function, among other things. 

We stress that the zeta function technology used on fractals applies equally 
well to Euclidean domains, manifolds, and graphs, and a fortiori produces the 
usual results on or U^. Furthermore, it unifies the treatments of massless and 
massive Bose gases: see 20 for a nice exposition in the W'- case. 



1.4. Main results. Below is a summary of the major results on the thermody- 
namics of massive and massless Bose gases in Sierpinski carpets. All terminology 
will be explained in subsequent sections. 

— Bose-Einstein condensation (Section [S]). For a low-temperature, high-density 
ideal massive Bose gas in an unbounded GSC F^o, Bose-Einstein condensation 
occurs at positive temperature if and only if the spectral dimension of the 
carpet ds{Fao) > 2. In this case, whenever the Bose gas density exceeds 
Pc(/^) := u^^sfkC (^): where Ci(i^oo) > 1 is a constant depending on Foo 



(see Proposition 5.3), then any excess density must condense in the lowest 
eigenfunction of the Laplacian. See Theorems |5.4| and |5.8[ 
Blackbody radiation (Section [6]). Let F be a GSC. At inverse temperature /3, 
the energy per unit spectral volume of photons inside LF is 

£{f3,L)^p-^''^^^^+^^HJ-log(^))+o{l) as L ^ 00, 



2L^ , 

where Hi is a periodic function of period ^diu{F) log I p given in Proposi- 



_ . . , 2 

tion 



6.3 and dw{F) is the walk dimension of the carpet. 
Casimir effect (Section [7|. Consider a waveguide Ha^b = o-F x [0, h] modelled 
after a 2-dimensional GSC F, and impose Dirichlet conditions on the outer 
boundary. If one places a pair of parallel plates at aF x {0} and aF x {6}, 
respectively, then the zero-temperature Casimir pressure on each plate is given 

by 

Pcas(a,5) = 6-('^=(^^)+i)ff2(^-logQ^^-Ho(l) as a ^ 00, 
where ds{fi) = ds{F) + l is the spectral dimension of the waveguide fig^b, and 



H2 is a periodic function of period ^dw{F) log If given in Proposition 7.2 
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These are to be compared with the classical "textbook" results in Euclidean 
space: 

— For a low-temperature, high-density ideal massive Bose gas in R'^ or Z'', Bose- 
Einstein condensation occurs at positive temperature if and only if c? > 3. In 
this case, whenever the Bose gas density exceeds Pc{(3) :— C (f); any 
excess density must condense in the lowest eigenfunction of the Laplacian. 

— At inverse temperature /3, the energy per unit Euclidean volume of photons 
in a cube [0, i]'* is 

^(/^' ^) = ^ ^(5^^ (^)c(rf +!) + «(!) asL-^oo. 

When d = 3 we recover the famihar "T" law," £{(3, L) = 7rV(30/34) + o(l). 

— Consider the rectangular waveguide [0,a]^ x [0,6] with Dirichlet boundary 
conditions, with [0,a]^ x {0} and [0, a]^ x {6} being the two parallel plates. 
Then the zero-temperature Casimir pressure on each plate is given by 

While our results are stated for true fractals, we believe that many of these 
behaviors can already be seen on finite-level approximations of the fractal, as 
indicated by our numerical work on the spectrum of the Laplacian (J. P. Chen 
and R.S. Strichartz, preprint [CS], M. Begue, T. Kalloniatis and R.S. Strichartz, 
arXiv:1201.5136 [BKS]). It would therefore be edifying if experimentalists can 
take on the challenge of constructing fractal-based quantum systems using, e.g. 
metamaterials, optical lattices, or superconducting qubits (for instance the pro- 
posal by [so]), and testing our results. 

The paper is organized as follows. In Section [2] we recapitulate several re- 
cently established results from the analysis on GSCs, namely the uniqueness of 
Brownian motion and the estimate of the heat kernel trace. In Section [3] we in- 
troduce the spectral zeta function on GSCs, show that it admits a meromorphic 
extension to C, and give its poles and residues. After recalling the rudiments of 
quantum statistical mechanics in Section|4j we then present our thermodynamic 
computations of the massive and massless Bose gases in Sections [5j |6] and [7] We 
conclude with Section [8] by discussing a special case of interacting Bose gas on 
Sierpinski carpet graphs, and offering some open problems. 



2. Established results on Sierpinski carpets 

2.1. Existence and uniqueness of Brownian motion. In this subsection we give a 
quick account of how the Laplacian on GSC is constructed. The precise details 
are high ly nontrivial and involve various techniques in potential theory (see 



e.g. 16 for background): we refer the reader to the original literature. For the 
purposes of this paper, it is enough to recognize the following mathematical 
facts. A Laplacian A is in 1-to-l correspondence with a nonnegative symmetric 
Markovian quadratic form, called the Dirichlet form £{u,v) = Ju{—Av), on 
an appropriate Banach space. Moreover, the said Laplacian generates a Markov 
process, which in our setting is either a simple random walk (on a graph) or a 
Brownian motion (on a subset of E**). 
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The Barlow-Bass construction [2j|5]. Let W^^ be a reflecting Brownian 
motion on the nth approximating domain of F. In order to produce a Brow- 
nian motion on the fractal F, one has to use the self-similarity of the carpet and 
take a suitable scaling limit. Barlow and Bass proved that there exists a family 
of time-scale factors {a„}„ satisfying 

<a^<c,^-j^) , (2.1) 

for some 01,02, Pf G (0, 00) independent of n, such that the sequence of sped- up 
Brownian motions X" = ^ on F„ has a subsequential limit. (The role of 
Pp will be discussed in the next subsection.) Any such limit process Xt, which 
respects the symmetry of F (henceforth referred to as _F-symmetric), is called 
a Brownian motion on F. We denote its infinitesimal generator by Cbb, the 
Barlow-Bass Laplacian (with Neumann conditions on doF). 

One can also rephrase the above result in terms of Dirichlet forms. Let 
Vn{dx) = {lp/mp)"dx be the Borel probability measure on Fn which assigns 
equal weight to each n-th level cell of F„. (Note that i^n converges weakly to 
a probability measure which is a constant multiple of the (i;i(i^)-dimensional 
Hausdorff measure, on F.) We introduce the Dirichlet energy on (F„, Vn) 

£n{u) = / \Vuix)f ly^idx) (2.2) 

■JF^ 

for all u e L'^{Fn,Vn) such that £n{u) < 00, and obtain the corresponding 
Dirichlet form by polarization: £„(u,w) = ^ [£n{u + v) — £n{u — v)]. Then the 
rescaled Dirichlet energies S'n{u) a„£'„(u) converge in subsequence to 

<S'bb(u) ^sup^{{l~Tt)u,u)L2/p^^), (2.3) 

where — e*'^^^ is the semigroup associated with Xf. The corresponding 
Dirichlet form is strongly local, regular, conservative, .F-symmetric, and self- 
similar: 

(^bb{u,v)= pfSbb{uo^s,vo^s)- (2.4) 

S£Si(Ft) 

The Kusuoka-Zhou construction |22 . Given i^„, let G„ = {VmEn) be 
the graph whose vertices lie at the centers of the cells 5" S iS„(F„), and whose 
edges connect vertices in nearest neighboring cells: see Fig. [2] Define the Dirichlet 
energy to be the usual graph energy 

^n{u) = ^ [u{x) - u{y)Y 

for all continuous functions u,w G C(G'„;M), and obtain the Dirichlet form via 
polarization. By checking a series of geometric conditions for which Poincare and 
Harnack inequalities hold, Kusuoka and Zhou were able to prove that {p^£n}n 
converges in subsequence to a strongly local, regular, conservative, i^-symmetric 
Dirichlet form S'kz satisfying the self-similar identity 

'?Kz{u,v)^ PF^Kz{uo<Fs,vo^s)- (2.5) 

seSiiFi) 
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The corresponding Brownian motion has infinitesimal generator CKZ^ which we 
call the Kusuoka-Zhou Laplacian. 

A question which lingered for almost two decades was whether the two con- 
structions yield the same limiting Laplacian on F . This was settled definitively 
by Barlow, Bass, Kumagai & Teplyaev. 

Theorem 2.1 Theorem 1.2]). Let F he a GSC equipped with the self-similar 
measure v. Up to scalar multiples, the set of non-zero, local, regular, conserva- 
tive, and F -symmetric Dirichlet forms on {F, v) contains at most one element. 



An equivalent statement to Theorem |2.1| is that up to deterministic time 
change, Lbb and Lkz both generate the unique F-symmetric Brownian motion. 
Henceforth we will denote this unique Laplacian by A. 

It should be noted that through extensive numerical computations [CS,BKS], 
we can demonstrate that both versions of the Laplacian on finite approximations 
Fn of F coincide at the bottom of the spectrum, as Fig. [3] shows. Notice that 
upon removing the power-law growth from the eigenvalue counting function, 
the remainder exhibits logarithmically periodic oscillation, which indicates the 
fractal nature of the geometry (see Section 2.3 for more details). For large enough 
n (typically > 3), the first [n — 2) segments of the two spectra, and to a 
lesser extent the (n — l)th segment, agree very well. Then deviation creeps in 
the nth segment: the Barlow-Bass spectrum becomes more characteristic of the 
spectrum on [0, 1]'', while the Kusuoka-Zhou spectrum is truncated due to the 
finite cardinality mj, of the graph. 

So despite our inability to state precise results rigorously on the finite approx- 
imations Fn, we believe strongly that all the asymptotic thermodynamic results 
presented in this paper, which are stated for F but invariably tied to the bot- 
tom of the spectrum, should be visible on F„ for sufficiently large n, beginning 
perhaps from n = 3. 



2.2. Resistance and heat kernel estimates. The aforementioned coefficient pp 
in (2.1 ) is called the resistance scale factor of the carpet F . When viewing each 
Fn as an electrical network, pp gives the renormalization factor relating the 
resistance of F„ to that of F„+i. To this date there is no known closed form 
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Weyl ratio of Barlow-Bass Laplacian on SC (a=0.87) 




- Level 4 

- Level 5 



Weyl ratio of Kusuoka-Ztiou Laplacian on SC (a=0.90) 




- Level 4 

- Level 5 




Fig. 3. Numerica l W eyl ratios W{s) = s~'^''^^N{s) associated with (a) the Barlow-Bass 
Laplacian [CS] and |(b)| the Kusuoka-Zhou Laplacian [BKS] on approximations F4 and F5 of the 
standard two-dimensional Sierpinski carpet F = SC'(3, 1), where N{s) := #{A £ Spec(— zi) : 
A/Ai < s} is the eigenvalue counting function of the Neumann Laplacian, normalized by the 
lowest nonzero eigenvalue Ai . 



expression for pp. The best known bound, obtained by cutting and shorting 
resistances, is [5| Proposition 5.1] 



mp 



^ <PF< 2^-'^If. 



(2.6) 



The connection between resistance and Brownian motion is as follows. Let 
dw{F) — \og{pFmp) / \oglp be the walk dimension of the carpet: this is the 
time-to-space scaling exponent for Brownian motion. For instance, if t{A) = 
mi{t > : Xt ^ A} is the exit time of a Brownian motion Xt from the set A, 
then dw is defined through E^[T{B{x,r))] x r''"', where is the expectation 
with respect to the law of Brownian motions started at x. 

A crucial point to make here is that while manifolds and Euclidean domains 
have dw = 2, fractals have d^ > 2, as can be seen through 6.17. (2.6). In a 



nutshell, this means that Brownian motion proceeds "more slowly" on fractals, 
due to the presence of obstacles at all length scales. One way to see this explicitly 
is through the heat kernel pt{x,y), which is the integral kernel of the Markov 
semigroup Tt = e*"^ , or cquivalently, the transition density of Brownian motion. 
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On fractals, regarded as a nietric measure space {F,d,v), the (short-time) heat 
kernel obeys sub-Gaussian rather than Gaussian bounds [4]: 



Pt(x,2/)xCt-'^''/'^-exp(^-c(^^^^^^y""j , te (0,1), (2.7) 

The leading power-law term is t"^''/*^™ , as opposed to the t~'^/'^ on d-dimensional 
Euclidean domains. This observation leads us to introduce the spectral dimen- 
sion of the carpet ds{F) :— 2dh{F)/dw{F) — 2logmp/ \og{ppmp). Note that 
for any GSC in 1 < ds{F) < dh{F) < d. 

The spectral dimension is a physically meaningful dimension because it is tied 
to the macroscopic behavior of Brownian motion. Recall that a Markov process 
on an unbounded state space is said to be recurrent if, with probability one, the 
process returns to the origin infinitely many times. Otherwise it is said to be 
transient, i.e., with positive probability the process leaves for infinity. Roughly 
speaking, on a homogeneous space, = 2 is the dimensional threshold above 
which Brownian motion is transient, and below which it is recurrent. Brownian 
motion at dg = 2 is typically also recurrent, but its behavior is more subtle than 
the ds < 2 case. 

Indeed, for Brownian motion in an unbounded GSC -Foo, Barlow and Bass 
showed that the process is transient (resp. recurrent) if ds{F) > 2 oi pp < 1 
(resp. ds{F) < 2 or pp > 1) [Sj Theorem 8.1]. The same dichotomy holds on the 
pre-carpet F [Hj Theorem 8.7]. This distinction has important consequences for 
the thermodynamics of the massive Bose gas, as we will explain later. 

To make comparisons, we mention in passing that all post-critically finite 



(pcf) fractals (for definition see 19 29]), e.g. any d-dimensional Sierpinski gasket 
with d > 2, have pp > 1 or d^ < 2, and support recurrent Brownian motion. 
Meanwhile, as an example of a random fractal, the incipient infinite cluster (IIC) 
at criticality of bond percolation on Z*^ with d > 19, or of spread-out percolation 



with d > 6, has ds = 4/3, which was proved recently by 21 . 



2.3. Sharp estimates of the heat kernel trace. In what follows we call K{t) :— 
Tr(e*^) = JpPt{x,x)di^{x) the heat kernel trace. Since the Laplacian on GSCs 
admits an eigenfunction expansion, one can write K{t) = X]j°=o , where 
{Xj} are the eigenvalues of —A. Equivalently, K{t) represents the probability 
that Brownian motion started at any point x returns to the said point x at time 
t. 



Using (2.7) we can already deduce that on GSC, K{t) x t-'^'^/'^- for smah 
t. This estimate is typically enough on simple spaces, but on deterministic self- 
similar fractals there is logarithmically periodic modulation on top of the power- 
law dependence in t, which is attributed to the discrete scale invariance of the 
space. Indeed, using renewal theorem type arguments, Hambly and Kajino sep- 
arately proved the following short-time asymptotics of the heat kernel trace. 

Theorem 2.2 ([l7)[T8l). Let F be a GSC, v he the self-similar measure on F, 
and A be the Laplacian on {F, ly) associated with either Dirichlet or Neumann 
condition on doF. Then 

K{t) = t-'^'^''^^ [G{- log t) + o{l)] asti 0, (2.8) 
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where G is \og{ppmp)-periodic, and is bounded away from and oo. 

Remark 2.3. Let N{s) := #{A < s : A e Spcc(— zi)} be the eigenvalue counting 
function (or the integrated density of states) of the Laplacian. Then K{t) = 
e~**diV(s), i.e., K{-) is the Laplace-Stieltjes transform of iV(-). Using the 
estimate K{t) x t^d.j2 ^^^^ easily obtains the Weyl asymptotics N{s) x gd.s/2 
for large s. 



In view of (2.8), it is tempting to go a step further and claim that 

N{s) ^ s'^"^^ [h(log s) + o(l)] as s^oD, (2.9) 
where h is log(/9i?mi?)-periodic an d is bounded away from and oo. There 



much numerical evidence that (2.91 holds on two-dimensional and thr ee- 
dimension al G SCs [7| CS, BKS], as Figure |3] demonstrates. Indeed, if (2.9) is 



true, then (2.8) follows immediately. However, as pointed out in '17', Section 4.2] 
and ^18^ Section 9] , the inverse Laplace-Stieltjes transform involves sophisticated 
Tauberian theorems which are very difficult to prove. This explains why we have 
elevated the role of the heat kernel trace over that of the (integrated) density of 
states. 

More recently, Kajino (in preparation) proved a sharper estimate of the heat 
kernel trace than ( |2.8[ ) by exploiting the full symmetry of the GSC and its 
boundaries. 

Theorem 2.4. Let F be a GSC, v he the self-similar measure on F, and A 
be the Laplacian on {F, v) with Dirichlet conditions on DqF. Then there exist 
continuous, log{ppmp) -periodic functions : M — > M for fc = 0, 1, • • • , d such 
that 

d 

if(t) ==^i-'*'=/'*™Gfe(-logt)+C'(cxp(-ci"a^r^)) as 1 10. (2.10) 

k=0 

Here dk := dh{F n {xi ^ ■ ■ ■ x^. = 0}) . Moreover Gq > and Gi < 0. 

Remark 2.5. It is believed that the same estimate holds on GSCs with Neumann 
conditions on doF , but the proof is more elusive. (N. Kajino, personal commu- 
nication) 

As the Gfe are periodic, we can expand them in Fourier series: Gk{x) = 
Epgz'5'fe,pe2'^P"/'°sfl^ where R ppmp = d^{F)loglF. We should note that 
very little is known about the functions Gk, except for the signs of Go and 
Gi. There is strong numerical evidence that Go is nonconstant, with (maxGo — 
minGo)/Go,o typically on the order of lO'^ for GSCs in or [cs]; but this 
has not been rigorously established. 



We point out three important features of the estimate (2.10) which are crucial 
to our thermodynamic computations in the sequel: 

1. The polynomial terms all have nonpositive exponents (— dfc/rfto)fc, where d^ 
is the Hausdorff dimension of the codimension-fc outer boundary of F. In 
particular, do — dh{F), di — dh(doF), dd-i = 1, and dd = 0. 
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2. Since the Gk are real-valued bounded functions, we have Gk.p = Gk,~p, and 

\Gk,p\ < ^ jl^°'^^\Gk{x)\dx < C < oo for ah p e Z. 

3. The remainder term decays exponentially in t, which is better than the power- 



law decay in the estimate (2.8). 

The first two features guarantee that the thermodynamic quantities (energy, 
pressure, etc.) associated with GSCs will be finite and meaningful. The last fea- 
ture, which is the most important of all, enables us to carry out the Casimir 
energy calculation unambiguously and without making any ad hoc regulariza- 
tion. 

Now recall that the heat kernel trace on Riemannian manifolds M C M'^ has 
the short-time asymptotics 



Going back to (2.10), it is easy to check that if _F is a GSC, then the heat 



kernel trace for LF {L > 0) is given by K{L-H) = if^(- \og{L'H))t-'^'/'^ + 



lower order terms as i J, 0. Comparing this against (2.11), we see that LF has 
an effective spectral volume of (72^'^'^ (47r)''=/^, where C2 ^ Go- In fact, we can 
use Cesaro averaging to define more precisely the spectral content of a GSC, in 



the same manner that the Minkowski content of a fractal string is defined 23 
This definition also appeared in fl]. 

Definition 2.6. Let F he a GSC. Then we define the spectral volume of LF 
(L > 0) to be 

V.{LF):=L''^{47rf'/^ lim — ^ Go{-\ogt)f = {Anf'/^GofiL"^ . 

rn-oo 71 log it Jfi-n t 

(2.12) 

We emphasize that the anomalous walk dimension d^; > 2 of the GSC results 
in its spectral volume being distinct from its Hausdorff volume. As will be made 
clear in later sections, all "extensive" thermal observables {e.g. energy, particle 
number) scale with the spectral volume, not the Hausdorff volume. 

Moreover we denote dk,p ■= 2 + ^^^^ • These are the complex dimensions 



of the GSC in the sense of M. Lapidus 23 , which reflect both the various spectral 



contents (via the real part of dk^p) and the discrete scale invariance of the fractal 



(via the imaginary part). Using this notation, (2.10) can be rewritten as 

d 

ii:(f) = ^^Gfe_p<-'^''.p/2-HC'(cxp(-ci"s;i^j) as UO. (2.13) 



3. Spectral zeta function on Sierpinski carpets 

The spectral zeta function of a self-adjoint Laplacian Z\ on a bounded domain 
is given by 
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This can be written as a Mellin transform of the heat kernel trace 



C/i(s,7) = 



r{s) 



t'e~^^K{t) 



dt 
T' 



(3.1) 



whenever t he ri ght-hand side is defined. In the case of a GSC with Dirichlet 
boundary, (2.101 impUes that the absicssa of convergence for C,a{-,i) is located 
at Re(s) = cLqxi — dh{F)/dw It is natural to ask whether this function can be 
extended to the entire complex plane, save for a countable (possibly infinite) 
number of poles. The following theorem, which was reported in (B. Steinhurst & 
A. Teplyaev, arXiv: 101 1.5485), serves as the linchpin for all subsequent results 
in this paper. 

Theorem 3.1. The spectral zeta function s i— >■ Czi(s,7) of the Laplacian on a 
GSC F , with Dirichlet conditions on doF, admits a meromorphic extension to 



Proof. Without loss of generality we suppose that (2.13) holds for t e (0,1), 
while there exist 03,04 > such that K{t) < 036^*^** for t > 1. Then 



a(s,7)r(s) = /i(s,7) + /2(s,7) +/3(s,7). 



where 



,/2 



dt 



fc=0 p£l 



l2{s,j)^j t'e-'"o(exp(^^ct-^ 
14 (s, 7)1 < C3 



dt 
T' 



t'e 

To compute Ii we use the identity 

r-l 



S„ lt^C4t 



dt 



t 



f\s^^,t^-pdt ^ r\.y^(~l)\ (-1)"7" 
Jn t Jn ^ n\ ' t ^n\{s + n- 



p) 



for Re(j3) > 0. Using linearity of the polynomial terms we get 



(-l)"7"Gfe,p 



n\(s + n — dk d/2) 

k=0pel,n=0 ^ '''P' ' 

In deriving this expression we interchanged the order of p-summation and in- 
tegration, which we justify as follows. For fixed 7 7^ 0, the summand indicates 
that the simple poles are (ife_p/2 — Ng. Away from the poles s H> /i(s,7) is 
holomorphic, since for each fc = 0, 1, • • • , d and each r g C, 
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is summable over p e N. When 7 = 0, we find 



Gk,p 



s - dk,p/2 



Next we let a := l/{dyj — 1), and make the bound 



r{s) 



< c 



1 



r{s) 







dt 
~t 



< Ce 



1" 



r{s) 



7(s,7) 



where 7(5, x) is the lower incomplete gamma function. Since 

_ 7™ 
W)^^"'^^ = ^ ns + m + l) 

is entire in both s and 7, l2{s,'y)/r{s) contributes an entire-in-s component to 
Ca(s,7). 

For I3, it suffices to know that for Re (7) > — C4, 



1-^3(5,7)1 < C3 



(7 + C4)-M - 
J7+C4 



C3|(7 + C4) T(s,7 + C4)| 



is finite, thanks to the holomorphicity of the upper incomplete gamma function 

s H- r{s, x) when a; 7^ 0. 

So finally we obtain, for Re(7) > — C4, the Mittag-Leffler decomposition 



^) = r(s) L L + „ _ 4 ^/2) 

This proves the theorem. □ 



+ (entire function in s). (3.2) 



Now we can enumerate all the possible poles and residues of the spectral zeta 
function on a GSC. 

Corollary 3.2. For a GSC F with Dirichlet condition on dgF: 

1. The poles of C,a{',j), where Re(7) > — C4, are contained in the set 



uuu 

k=OpeZ n=0 



dk,p 



with corresponding residues 

Res ( CA{-,j),-n + 



dk,p 



(-l)"7"Gfc,p 
n\r{-n + dk,p/2y 
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(a) 
Im(s) 



(b) 
Im(s) 



Re(s) 



1 1 
2 





• • 




• • 


d^,(F) log Ip 


• • 





_L 








• • 

• • 

• • 



Re{s) 



Fig. 4. The singularities of s i-^_Ca{s,0) associated to (a) a 2-dimens iona l parallelpiped, 
using the Epstein zeta function; |(b)| a 2-dimensional GSC, usmg Corollary |3.2| A hollow circle 
indicates that the singularity is either a removable singularity, or a pole with zero residue. 



2. The poles o/C^(-, 0) are contained in the set (J 
ing residues 



, with correspond- 



Res (^a(-,0), 
)) is ar 

3. Czi(s,0) = for all s e -N 



2 J r{dk,p/2)- 

In particular, C/l(5,0) is analytic for Re(s) < 0. 



Proof. The poles and residues can be read off from (3.2), so the first two claims 
are immediate. Note that dd.o/^ — N = — N are poles of Czi(', 7) with zero residue 
because l/r{s) — for all s £ — N. The last claim follows from the observation 
that when Re(s) < 0, (a{s, 0) is a product of l/r{s) and a holomorphic function. 
□ 

Fig. [4] depicts the singularities of Czi('jO) for a typical GSC versus that for 
a parallelpiped. In both domains, the poles with the largest real part represent 
the spectral volume dimension. What distinguishes a GSC from an Euclidean 
domain is the "tower" of poles, or complex dimensions, along the imaginary 
direction at each Re(s) = dk,o, with even spacing of 27ri/((i^ log /i?). Thus the 
discrete scale invariance of the GSC manifests itself in two ways: it produces 
log-periodic modulations of the heat kernel trace and, via a Mellin transform, 
creates towers of complex-valued poles in the spectral zeta function. 



4. Reviev^r of quantum statistical mechanics 

Here we briefly review notions of quantum statistical mechanics that are needed 
for the rest of the paper. Much of this is known to physicists, but we use the lan- 
guage of functional analysis to make notations precise. For a complete treatment 
please consult, e.g., [sj Chapter 5]. 
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4-.1. Gibbs state and partition function. The Bose (resp. Fermi) gas is a system 
of many quantum particles satisfying an appropriate rule of quantum statistics. 
Let us denote the Hilbert space for a single quantum particle confined to domain 

C M"^ by Hi := L'^{F): the state of the particle is described by a wavefunction 
ip & T-Li- For a system of n such identical bosons (resp. fermions), the wave- 
function ^n{xi, ■ ■ ■ TXn) € L^{F'^) is (anti)symmetric under particle exchange: 
V'nC' • • - ■ ■ ,Xj,- ■ ■) = ±'0„(- • • , Xj, ■ ■ ■ ,Xi,- ■ ■). Therefore the appropriate 
n-body Hilbert space, denoted by T-Ln, is the (anti)symmetric subspace of Hf". 
By default we set "Ho = C. 

To describe the equilibrium thermodynamics of the quantum many-body sys- 
tem, it is convenient to adopt the grand canonical ensemble, where the particle 
number is not kept fixed. Mathematically, we take the Fock space T = Hn 
as the underlying Hilbert space. For each one-body self-adjoint operator W : 
Hi^ni, define Wn'-Hn^ Tin by 

{0, if n 

n 
^P±(^i®...0(W^^,)^---^V'n),ifn>l 

for all ijjj G Hi, where P± is the (anti)symmetric permutation on the components 

of the tensor product. Then we call dr{W) := second quantization 

of W on T. An important example is where W is the identity operator 1 on 'Hi: 
then 1„ is n times the identity operator on Hn, so its second quantization reads 
dr{l) = ©^o"" ~- ^^^'^ known as the number operator on T. 

For our purposes it suffices to consider a one-body Hamiltonian operator 
H which is bounded below. Let (3 e (0, oo] be the inverse temperature, fj, € 
(— oo, inf Spec(iy)] be the chemical potential, and H := dr{H). Consequently 
:= e^^'^'"'^^^^^^ is a trace-class operator on J^, and its trace 

^p,, Tr^P^,. = Tr^e-«H-^^) 

is called the grand canonical partition function of a thermal system with Hamil- 
tonian H . The free energy of the system is 

The Gibbs state i^p.fi at {13, ii) is a linear functional over the C*-algebra on T 
such that 

Physically, if A belongs to the (*-)algebra of observables on Hi, and A = dr{A), 
then w^,^(A) represents the expectation value of A in the grand canonical Gibbs 
state. To experts of operator algebra, we remark that the Gibbs state is a (r, /?)- 
Kubo-Martin-Schwinger (KMS) state [sj Definition 5.3.1]: namely, for all A,B 
in the algebra of observables on J^, 

c.^,^(Ar,0(B)) = o.^,^(BA), where ^(O) :- e'*(H-^N)og-.t(H-MN)^ 

Two observables we will be interested in the sequel are the expected number 
of bosons _ 

^pA^) = r^^^ = logs,.,, 
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and the expected energy 



1 5^/3,/. d „ 



For a non-interacting Bose gas the expression for the partition ftmction can be 
reduced to a trace over the one-body Hilbert space Hi- Assume that H has pure 
point spectrum with increasing eigenvalues < Eq < Ei < ■ ■ ■ < Ei < ■ ■ ■ ^ oo. 
Then 

i ni=0 i 



Therefore 



logS^,^ = -Tr„, log(l - e-^(^-'^)), 



The following fact is useful and will be invoked several times: /i h- > cj^^(N) is a 
convex, monotone increasing function which maps (— oo, inf Spec(if)] to (0, oo]. 



4-2. Thermodynamic limit. Since we will deal with the delicate topic of phase 
transitions on fractals, it is important to discuss how we take the thermody- 
namic limit. As is customary, we exhaust an unbounded space Soc by an in- 
creasing family of finite subsets {S'„}„ with suitable boundary conditions. The 
thermodynamic limit of an "intensive" thermal observable {e.g. pressure, par- 
ticle density) on Sco is then understood as the limit of a sequence of the said 
observables on Sn, subject to proper constraints. 

For the unbounded carpet Foo, we choose the natural exhaustion := 
{lpF}n. This involves scaling up the size of the system by Ip each time, which 
in effect scales down the Laplacian by Ip^. To make this length dependence 
explicit, we will denote by A the dimensionful Laplacian on a bounded domain 
not necessarily of unit characteristic length. If the characteristic length is L, then 
we write zi^, which is L^^ times the dimensionless Laplacian A (on the domain 
of unit length). In the same manner we will add a subscript L to any quantity 
which depends explicitly on the domain size L. 

We will show shortly that, due to the log-periodic modulation inherited from 
the heat kernel trace, the thermodynamic limit along the natural exhaustion 
does not exist in the literal sense, but rather in the weaker limsup/liminf sense. 
Of course our choice of exhaustion is by no means unique. There might be other 
choices of exhaustion for which a stronger limit is attained, btit we suspect that 
those would marginally improve the results to be presented below. 
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5. Massive Bose gas in Sierpinski carpets 

In this section we focus on the Bose gas with H = — A and fi < inf Spec(iJ). This 
Hamiltonian captures the behavior of a gas of non-relativistic, non-interacting, 
and charge-neutral atoms. The partition function of the massive Bose gas in a 
general bounded domain is given by 

log Sp,^ = -Trn, log (l - e-''^"^-^)) . (5.1) 

Its connection to the spectral zeta function is the following. 

Proposition 5.1. The grand canonical partition function S"^,^ for an ideal mas- 
sive Bose gas in a bounded domain F is given by 

1 r'^+i°° d 
\ogSp^^ = — r'r{t)at + l)CA{t,~^i)dt, a>^. (5.2) 



Proof. Using the T aylor expansion for log(l — a:) and functional calculus, we 
expand the RHS of (5.1 1 in a power series: 



CO _^ 



-n0{-A-fi) 



n=l 



By the identity 



we have 



1 
27ri 



a~*r{t)dt, cr > 0, 



1 ^ 1 /"CT-i-lOO 

logS^.^ = Tr„,— ^- / {nl3r'rit){-A-fi)-*dt. 

It remains to interchange the order of integration and summations. To do so 
we need to pick a such that both Y.'^^i n-^''+"'+^^ and Tr^^(-^ - 
converge for any c € M. This happens when a > W ^ = □ 



5.1. Bose-Einstein condensation in the unbounded carpet. For the massive Bose 
gas we are interested in the density of bosons, that is, the number of bosons per 
unit spectral volume: 



c.L,ff,^(N) 1 1 1 d ^ 

^^(^'^^ V.{L) - K(Z) ''^e/'C-^.-A.)-! - WMd'^ ^^^'-'^''' 

(5.3) 

Lemma 5.2. Let F be a GSC. Then the density of ideal massive Bose gas in 
LF at {13, /i) is 

(47r) 2 Go,o \fc=Opez/ P m=o V ^ 



(5.4) 

where (X]fc=oSpez) "means the double sum excluding the term {k,p) — (d, 0) 
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Proof. W e ev aluate the contour integral in (5.2 1 via a residue calculation. By 
Corollary : 



Upon exc 



3.2 



Czv(-,7) has poles {jnenaidk.p/^-n) with residues Jr,!,^^^ 



uding the poles with zero residue, — No, one finds 



(-l)"7"Gfc.p 
(-n+dfc,p/2) • 




The lemma follows from taking the /i-derivative of log Sl,i3,h and then dividing 
by the spectral volume Vs{L) = (At:)^ Gq^qL'^' . □ 

We are in a position now to take the thermodynamic limit of the particle 
density along An -Foo- To streamline the arguments, we introduce the fugacity 
z = e'^^, and state all quantities in terms of z rather than /i. We write pA^iP, z) 
for the particle density in An at (/3, z). 

Proposition 5.3. Let 'pdP) = limsuppyi„(/3, 1) and p (/3) ~ liminf /5yi^(/3, 1) 

TL — ^oo n—^oo 

he, respectively, the upper and lower critical density of the massive Bose gas in 
Foo. Then pj^(3) = = oo i/4 < 2, while < pjy(i) < p^(/3) < oo i/d, > 2. 

In the latter case, 

- iR\ 1 maxGo /ds\ , n-. 1 minGo /4\ ,^ 

(47r/3)— Go,o V ^ / (47r/3)— Go,o V ^ / 



Proof. First we note that as L — > oo, only the volume terms (fc = 0) dominate: 

"'•A^' ^ E «... £ 

(47r/3) 2 Go,o pez ^ ^ m=o V ^ / 

By expanding the Riemann zeta function in a series and using some manipula- 
tion, we arrive at a more transparent expression 

(47r/3) 2 Go,o ™=i V \^ J J 

Replacing L with Vp we conclude that as n — > oo, 

(4^/3)tGo,o ^1 V \{lFr'JJ 

from which the proposition follows. Finiteness comes from the convergence of 
the Riemann zeta function C,{ds/2) = X]m=i "1"^°^^. Note that we do not claim 
that p {(3) — PciP) since Go is probably nonconstant. □ 
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The consequence of a finite upper critical density is that any excess Bose gas 
must occupy the lowest eigenfunction of the Hamiltonian, a phenomenon known 
as Bose- Einstein condensation (BEC). This is the content of the next theorem. 
We will state the result for a fixed-density Bose gas, i.e., fix the particle density 
Ptot along An Foo- This of course requires adjusting the fugacity Zn for each 
n. Moreover, we denote by Ek{L) and Ek{An) the (fc + l)th eigenvalue of the 
Laplacian on, respectively, LF and An. 

Theorem 5.4. Assume ds > 2. For each ptot > 0, let Zn be the unique root of 

PA„W,Zn) = Ptot- 

(i) If Plot ^ ~Pc{P)j cindz is the root of limswp py[^{/3,z) = ptot, then liminf 2:„ = 

n — ^oo 

Z. 

(ii) If Ptot > T>c{P)) then lim z„ — 1. Moreover if we denote the boson occupation 
density in the ground state, or the condensate density, by 

then lim [p^ (/3, z„) + p^„ (/3, 1)1 =ptot- 

The proof, which mirrors that of [sj Theorem 5.2.30] with some modifications, 
relies upon extensive use of the convexity of z i— >■ pl(/3, z) and the lemma below. 

Lemma 5.5. For any zi > Z2, 

Pl{P,Z2) Pl{I3,Zi) ~ Pl{/3,Z2) PlW^Zi) 

Z2 - ZI-Z2 - ziil-zie-P'^o{L)y ^^'^^ 

Analogously, if we denote the boson occupation density above the {m + l)th eigen- 
function by 

k>m ^ 

then for any zi > Z2, 

pT{I3,Z2) ^ Pt+(j3,z,)~p^[:+{f3,z2) Pr(A^i) o^ 

Z2 - Z1-Z2 - zi(l-zie-/5^™W)- ^ ' 



Proof. Using the convexity of z i— >■ pl{P, z), we have for any zi > Z2, 



^P^tR ^ ^ Pl{I3,zi) - pl{(3,Z2) dpL , . 

-^(P,Z2) S S -^[P,zi), 

OZ Zi — Z2 oz 



where the derivative is easily computed: 
dpL . 1 



dz 



.z) 



Vs{L) 
1 



-Tr 



Tr?ii 



1 



-lg-/3/i^ - 1 1 - zc^Al 



(5.9) 



(5.10) 
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Since —Aj^ > Eq{L) and z < e ^—^ by assumption, we can make the bound 



< 



1 _ ze^^i^ ~ 1 - ze-^^o(L) ■ 



(5.11) 



Applying (5.10) and (5.11) to ( |5.9| ) yields (5.6). The same argument carries over 
to p™^ by changing Eo{L} to Em{L) in the upper bound. □ 



Proof (of Theorem 5.4)- Since Go > and Gi < 0, we have for all sufficiently 
large n, 



1 



< 



(4^/?)^Go,o 
1 

(47r/3)^Go,o 



£z"Go(-log 



m— 1 

oo 



m/3 



m 2 



z^m"^ max(Go) = limsuppyi„(/3, 2;). 



Now we prove (i) By the preceding inequality (5.12) and the convexity of 
Pl(/3, ■), we deduce th at Z n > z for a ll su fficiently large n. Then together with 
the first inequality in (5.6 ) of Lemma 5.5 with zi = z„ and Z2 = z, we get 



< Zn — z < z 



PA„W,Zn) 

PA„{l3,z) 



-1 =z 



Ptot 



- 1 



Taking the liminf of this inequality yields the result. 

To prove (ii) , suppose ptot > 'pAP)- If -^n < 1 for some large n, then Pj,(/3) < 
pyi^(/3,z„) < limsupp/i„(/3, z„) < limsuppA„(/3, 1) —'■ PciP)i which is a contra- 

n— >-oo n— >oo 

diction. Thus z„ > 1. But since z„ < e^-^°'^^"^ we have lim z„ = 1. Moreover 

n— ^00 

p™+(^, zl) > p™"'"(/3, 1) for any rn € No, so in particular 



p°+(/3,z„) >p^+(/3,l) >p^„(/3,l) 



^0+ 



1 



VAAn)l3EoiAn)'' 



and 



z„) + p^„ 1) < p^„ z„) + p°+ z„) + 



1 



VMpEoiAr^ 



Ptot 



Noting that [y3(A0^^o(^«)]"' = Gn^-'^^ 



as n — >■ 00, we ffiid 
liminf [pV (/?, ^n) + PA„(/3, 1)1 < Ptot- 

n— foo 

It remains to show that lim sup [p/i„(/3j ^n) + PA„(/3, 1)] ^ Ptot- Here we uti- 

n— f 00 

lize the fact that for any m G N and n E Nq, the occupation density in the 
(m + l)th eigenfunction of the Laplacian on An 



1 



1 



< 



1 



1 



Vs{An) Zn^ePE^i^-) _ 1 - Vs{An) eMEm{A^)-Eo{A^)) - 1 

1 1 
- VMn) P{E,n{An) - Ea{An)y 
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The right-hand side scales with '^'^ and hence tends to as n — > cx), indepen- 
dently of TO. This implies that 

Ptot= lim [p5iJ/3,z„)+p°+(Az„)] = lim [pO J/?, z„) -t- ^n)] 



for any m G N. By replacing L with yl„ in the second inequality in (5.8) of 
Lemma |5.5[ and setting zi — z„ and Z2 = 1, we obtain 



It follows that for all n and for all sufficiently large m, 



where the Weyl asymptotics of the Laplacian was used in the last line. Thus 



Ptot < lim sup 



p°J/3,z„)-Hp^„(/3,l) (i-2Cto1-^)" 



As the left-hand side is independent of to, we can take the limit to — >■ cx) to 
complete the proof. □ 

Remark 5.6. Here we encounter a subtlety not seen in the Euclidean setting, 
namely, that the length scale enters into the log-periodic modulation Gq in the 
particle density. Therefore the thermodynamic limit needs to be handled with 
care. Theorem |5.4| says that in the lim sup density sense, condensation occurs 
if Ptot exceeds the upper critical dens ity p ^(/3). This appears to be the optimal 
threshold as our next result (Theorem |5.7[ ) suggests. Notice also that neither the 
condensate density p\ (/3, z„) nor the critical density p^^ (/3, 1) likely has a limit 
as n — > oo, but the sum of the two converges to ptot- 

One might ask whether a stronger sense of convergence could be used to char- 
acterize the thermodynamic limit on fractals. A possibility would be to use the 
Cesaro averaged density in place of the lim sup density. However incorporating 
the former into the convexity proof of Theorem |5.4| appears difficult. 



We have now shown that for an ideal massive Bose gas in an unbounded GSC 
with da > 2, there is macroscopic occupation of particles in the ground state at 
sufficiently low temperatures. It is well known that Bose-Einstein condensation 
is a quantum phase transition. To illustrate this point, we will show that in 
the thermodynamic limit L — > oo, the free energy density of the Bose gas 

becomes non-analytic in an appropriate sense. This is a familiar criterion of phase 
transition in the physics literature. A more technical criterion, namely, the non- 
uniqueness of Gibbs (or KMS) states at low temperatures, will be discussed 
elsewhere. 
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A straightforward computation shows that as L ^ oo 
1 



(4^)-/3^+iGo,o ^1 



log 



771/3 



m-C^+i) +o(l). 



Taking the sequence of free energy densities along yields the following 

result. 

Theorem 5.7. For each ptot > 0, let Zn be the unique root of pA„{(3, Zn) — ptot- 
(i) If ptot < Pc{l^)j andz is the root o/limsup/3/i„(/3, z) = ptot; then 



lim inf /^^ Zn) = lim inf fA„ iP, z) 



(47r)*/3T+i Go 



max Go . 



(5.12) 



where Li^ (z) — > — zs i/ie polylogarithm. 



r^y) //ptot > Pc(/5). <^en lim [/^„(/3,z„) - /.i„(/3, 1)] = 0. 

n— *-oo 

Proof. The key is to realize that ~-fL{l3,z) > 0, and that z i-^ — z) is 
increasing and convex. From this it follows that for any zi > Z2, 



fL{l3,Z2) - fL{l3,Zi) dfL,„ . 

Zl — Z2 OZ 



1 



Pl{I3,zi). 



(5.13) 



For (i) we replace L, zi, and Z2 in (5.13) with An, Zn, and z respectively to 
find 



0</.i„(^,z)-/^„(^,z„) < 



Ptot. 

/3 



1 - 



Take the lim in f on all sides and use Part |(i)| of Theo rem 5.4 to finish. 

Similarly for (ii) we replace L, zi, and Z2 in (5.13) with A„, Zn, and 1 respec- 
tively. By virtue of Part |(ii)] of Theorem 5.4 we can take the limit this time to 
deduce the result. □ 



Part |(ii)| of the preceding theorem says that in the condensation regime the 
thermodynamic free energy density is, in a proper sense, independent of the 
particle density ptot, or of the (nonzero) condensate density. This is because in 
the thermodynamic limit, the condensate resides in the zero-energy eigenfunction 
and does not contribute to the free energy. For example, in the sense of Cesaro 
averaging, we have 



n 1 ™ 

lim -^/^J/3,z,)= lim -Y^f^^ip,!)^- 



1 



n— >oo ri 



fc=l fc=l 

On the other hand, in the liminf sense 



1 



(4^)^/3^+1 
max Go 



1 



liminf /^„(/3,z„) = liminf /a„ (/3, 1) = — d,^, — 5, Lid^ (1). 

n^oo n^oo (47r)^/3^ + l Go,0 ^ 



Combined with the result in Part 



(i) one may check that liminf /yi^ (/3, z„). 



regarded as a function of ptot G (0, oo), is non-analytic at pJyjS), since the left 
derivative (a nonzero constant times C (^)) does not equal the right derivative 
(0). 
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MS(3, 1) M5(4,2) 


MS(5,3) 


M5(6.4) 




logg 20 Si 2.73 log4 32 = 2.5 


logs 44 ^ 2.35 


logg 56 2.25 


Rigorous bounds on 

ds 5; 


2.21 ~ 2.60 2.00 ~ 2.26 


1.89 ~ 2.07 


1.82 ~ 1.95 


Numerical ds [CS] 


2.51... 


2.01... 




BEC exists? 


Yes Yes 


Yes (?) 


No 



Table 1. A sequence of Menger sponges whose spectral dimensions c ross over 2. The criterion 
for existence of BEC in an ideal Bose gas is determined by Theorem |5.8| 



5.2. Connection to Brownian motion. 

Theorem 5.8. Let F he a GSC. Then the following are equivalent in the un- 
bounded carpet Foo •' 

1. The spectral dimension dg > 2. 

2. Brownian motion is transient. 

3. Bose-Einstein condensation occurs for a low-temperature, high-density ideal 
Bose gas. 

Proof. For Fqo, the equivalence between (1) and (2) is give n in [5) Theor em 8.11. 



The equivalence between (1) and (3) is implied by Lemma |5.2| and Theorem 5.4 
□ 

On general graphs, the connection between spectral dimension and BEC has 



been mentioned in previous literature [9 10 . The equivalence between tran- 
sient Brownian motion and BEC, on the other hand, appears to be a more 
general and robust criterion. To the author's knowledge, outside of and Z'', 



an explicit statement about the latter has only been made by 13 14 25] in the 
context of inhomogeneous graphs. In these works, the authors considered the 
adjacency operator, as opposed to the Laplacian, on various exotic graphs [e.g. 
comb graphs, star graphs) and density-zero perturbations of Cayley trees. (We 



will discuss their model briefly in Section 8.1 ) They gave conditions for BEC 
by analyzing the spectral properties of the adjacency operator, in particular its 
Perron- Frobenius eigenvector, and showed that the existence of "hidden spec- 
trum" implies the finiteness of the critical density. These developments lead us 
to believe that the deep connection between Brownian motion and BEC holds on 
very general spaces, including graphs, manifolds, and fractals. A careful study 
of KMS states on these spaces may uncover more details. 

We end this section by giving specific examples where one may probe BEC 
in non-integer dimensions. Recall that if is a GSC embedded in W^, then 1 < 
ds{F) < dh(F) < d. Also, an unbounded GSC whose cross-section contains a full 
copy of supports transient Brownian motion, i.e., the carpet has dg > 2 5 . 
So if one's goal is to find a GSC whose spectral dimension is close to 2, fne 
most natural candidates are the Menger sponges in M'^. For example, consider 
the family of Menger sponges {MS{a,a — 2)}a>3, the first four members of 
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which are shown in Table [1] AU such sponges have dh > 2, but only Af5(3, 1), 
MS'(4, 2), and possibly AIS{5,3) have dg > 2, in which condensation of ideal 
Bose gas is possible. In particular MS{5,3) lies just a hair above the critical 
dimension 2, as our numerics show. It would be very interesting if one could find 
a one-parameter family of fractal sponges whose spectral dimensions vary across 
2, and study the onset of BEC in such spaces. 



6. Blackbody radiation in Sierpinski carpets 

In the next two sections we discuss the Bose gas with H = \J —A and ^ — 0. This 
Hamiltonian corresponds to that of a gas of massless, non-interacting, relativistic 
particles, the prime example being photons. Here we have set /i = because 
the particle number cannot be kept finite in realistic experiments. The grand 
canonical partition function S"^ in this case reads 

logS^, = -Tr^,,log(l-e^''^), (6.1) 

where I) — \f--A_ is the dimensionful Dirac operator, carrying units of inverse 
length. 

Proposition 6.1. The grand canonical partition function Sp for an ideal photon 
gas in a bounded domain F is given by 

\ogSp = ^ / r'*/^(2i)C(2t + l)a(i,0)di, a>^. (6.2) 



The proof is essentially identical to that of Proposition |5.1[ so we omit it. 

Lemma 6.2. Let F be a GSC. Then the grand canonical partition function Sp 
for an ideal photon gas confined in LF is given by 

iog^.,. = 2|^^l (i) 7^;g^c(4. + i)G..+fa.(4o), 

(6.3) 

where ^X]fc=o X^pez) rneans the double sum excluding the term {k,p) = (d, 0). 




Proof. We evaluate (6.2) via a residue calculation. The first term of (6.3 



comes 
because 



from the residues at the poles of Ca(^'O)' though we exclude dd,o = 
it has zero residue (as well as being a simple pole of r{2t) and of (^{2t + 1) 
simultaneously) . Then there are contributions coming from the nonzero poles of 
r{2t), which are -N/2. But CA{t,Q) = when t e -N by Corollary while 



C(2< + 1) = when i G - (N +h)- Only the residue at t = -1/2 contributes and 
produces the second term of ( |6.3| . □ 

The first term of ( |6.3[ ), whose summands are all proportional to some non- 
negative power of L//^, represent thermal (positive-temperature) contributions 
from the various codimension-fc spectral volumes. On the other hand, the second 
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term, which is Unearly proportional to /3/i, has no thermal origins. To tackle 
this issue, let us add to (6.3) 



\ogS 



Cas 



(6.4) 



so we retain only the thermal contributions. Notice that this effectively adds a 
(free) energy 

Ec.s = -/3-MogHcas = \c,a(-\,0\ . (6.5) 



to the existing vacuum energy Hq, which was by default. (6.5) is the renowned 



Casimir energy representing the energy of vacuum fluctuations. It is half the 
regularized sum of the eigenvalues of the Dirac operator, where regularization 
should be understood as the analytic continuation of C/i('SjO) to s = —1/2. By 
Theorem |3.1[ such continuation is always possible on GSCs. The Casimir energy 
has measurable consequences and will be the focus of the next section. 

The most important phenomenon associated with thermal photon gas is black- 
body radiation. Here we are interested in the energy per spectral volume 



f(L,/?) 



1 



d 



Vs{L) 

as well as the isotropic pressure 

d 



Vs{L) dp 



log. 



'■L.,p,0, 



P{L,P) 



1 d 



Proposition 6.3. Let F be a GSC. As L — > 00, the mean energy density and 
isotropic pressure of a thermal photon gas inside LF are 



FiL,P) = 



oil), 



(6.6) 
(6.7) 



where Hi{-) is ^ log R-periodic with 



Hiix) 



ds ^ / d.s 



dg + l 

TT 2 



(ids + 1) 



1 



TT 2 Cxo,0 pgZ\{0} 



do.; 



C{do,p + l)rfo,pG'o,p 



(6.8) 



Proof. Straightforward differentiations of the thermal terms in (6.3), combined 
with the gamma function identity 



Fjz) ^ r{{z + l)/2) 
r{z/2) 2i-^0r ' 
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Fig. 5. The Sierpinski carpet waveguide a{GSC) X [0,b]. 



yield the results. We note that for aU s e M, \r{s + it) \ ~ V27r|t|''-i/2g-7r|t|/2 ^ 
t — >■ ±00 by Stirling's formula, while for s > 1 and t g M we have |C(s+zt)| < ^{s). 
These observations, plus the boundedness of the |Go.p|, are sufficient to ensure 



the summability of (6.8). □ 



We have intentionally split Hi in order to make comparisons with the Eu- 



clidean case. If only the first term in (6.8 1 was considered, then (6.6) would be a 



direct extension of the Stefan-Boltzmann law from Euclidean to noninteger di- 
mensions. But the discrete scale invariance of the GSC gives rise to log-periodic 



modulations depending on the temperature as well, viz. the second term in (6.8). 
Note also that the ^ in the periodicity ^ log R of Hi comes from the Dirac op- 
erator being the square root of the Laplacian. 

Meanwhile, ( |6.7[ ) represents the fractal version of the equation of state for the 
photon gas T. This is to be compared with the Euclidean version P = E/d, 
where d is the Euclidean dimension (also the spectral dimension, since d^ = 2). 



7. Casimir effect in Sierpinski carpet waveguide 

In this section we derive the radiation pressure on some boundary faces of a 
compact fractal domain, due to either vacuum (zero-temperature, j3 — 00) or 
thermal (nonzero-temperature, j3 < 00) effects. 



7.1. The parallel-plate setup. The original calculation, due to H. Casimir in 1948, 
concerned the vacuum energy in [0, a]^ x [0, b] for fixed a ^ 1, and the resultant 
pressure on the two parallel faces [0, a]^ x {0} and [0,a]^ x {h] as a function 
of b. The conventional method is to carry out a sum over the (exactly known) 
eigenvalues of the Dirac operator, take the a ^ 00 limit, then apply "zeta- 
regularization" to the formally divergent sum to produce a finite value. Similar 
calculations have also been done on other Euclidean domains or manifolds, such 
as spheres and cylinders. 

In the spirit of the classical parallel-plate setup, we will work with a Sierpinski 
carpet waveguide fla.b a{GSC) x [0, 6], where we take GSC c to be some 
2-dimensional Sierpinski carpet and fix a 3> 1. See Figure[5] The parallel plates in 
this case correspond to a{GSC) x {0} and a{GSC) x {6}. For consistency, we will 
impose Dirichlet conditions on the outer boundaries 9([0,a]^) x [0,5], a{GSC) x 
{0} and a{GSG) x {b}, and Neumann conditions on the inner surfaces. The 
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Hausdorff and spectral dimensions of the waveguide are, respectively, dh{f2) 
dh{GSC) + 1 e (2,3) and d,{(2) = ds{GSC) + 1 e (2,4(r2)). 



7.2. Casimir pressure at zero temperature. Our main objective is to compute 
the spectral zeta function Cn (^^jO)) or twice the Casimir energy. Because the 
Laplacian on Q is the direct sum of the Laplacian on GSC and the Laplacian 
on [0,6], the associated heat kernel trace factorizes: 



Kn^^^{t) = A'a(GSC)(i)A>,fc](t) = Kasc (a^'i) ^[oa] ip-H) . 
Recalling that the spectrum of the Dirichlct Laplacian on the unit interval is 



U^i {j^TT^}, we find 



OO ^ 



— -- + — V 



.-fit 



/Ant 



Vnt 

.1=1 



where the Poisson summation formula was employed in the third equality. Then 
using (2.101 with = 2 for Kcsc{t), we get 



EE^..(;;,V'%(:^fexpf-cfi 

k=0 pSZ 



/Ant 



(7.1) 



The spectral zeta function for the Laplacian on fia,b is the meromorphic exten- 
sion of 



n 



^ f°° df 

— / t^K,,^^{t)e'^*-, Re{s)> 
(s) Jo t 



ds{f2) 



(7.2) 



Theorem 7.1. As a — > oo, 



4-« 



1 V^Y^A /a\«'=f 



fc=o pel 



'^■"6-C(2 + d.,p)r(l+'f 



o(l) 



(7.3) 



Proof. Though the proof is similar to that of Theorem 3.1 we will be more 
precise in the computation here. Let /(<) = O ^exp (^—Ct~ By defi- 

nition, there exist T, Af e (0,cx)) such that \f{t)\ < Af exp (^-Ct~3^r=T^ for 
ah t e (0,T). It follows that \f{a-'^t)\ < M cxp (-Ca^ii^t^^^) for all 



Statistical mechanics of Bose gas in Sierpinski carpets 



29 



t € {0,a^T). So when computing the Melhn integral we spht the domain of 
integration (0,oo) into two intervals, (0,a^T) and [a^T, oo). 

Expanding the product in (7.11 and then plugging it into (7.2), we find 



1 



1 



I{t)e 



-It 







dt 
~t 



I{t)e 



dt 

T' 



where I{t) contains terms of the following types: 



1. Polynomial terms of the form t^^, Re(p) > 0. 

2. Product of a polynomial term and an exponential term: exp(— csi^''), 
Re{p),r] > 0. 

3. Product of a polynomial term and two exponential terms: 
t-Pexp(-C6t"'') exp(-C7i"i), Re{p),ri > 0. 

For the polynomial terms we have 



dt 

^-Pg-7t_ - ^,P^.f 



t 



where 'y{s,x) and r{s,x) are, respectively, the lower and upper incomplete 
gamma functions. Then we need to use the asymptotics of the incomplete gamma 
functions near the branch point x = 0: x~'^-^{s,x) — >■ 1/s and, for Re(s) < 0, 
x^^r{s,x) — l/s. These imply that 



lim 

7->-0 



dt 

t-Pe-^'- 



°° dt 



t 



= 0. 



Next we tackle terms of type (2). By making a change of variables we find 



Ji(a,6,7) := 



dt 



/,! du 
u 



Notice that for u G {c^ / {a^Ty ,oo) we have the upper bound 
gl7k g^j^jj ^jjg remaining integrand is integrable: 



-7(c5/n)i/'' 



du 



r 



/c5/(a2T)') ' 

So by the dominated convergence theorem 



P C5 



< oo. 



lim Ji(a, 6, 7) 
7— >o 



1 



The same argument applies to 

J2(a,6,7) := 



yP/')g-« 



du 



-r 



p C5 



2rp 



t-Pe-'''"'e-^*j, 



< 
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yielding 



Thus 



lim 

7->-0 



lim J2 (a, b, 7) = ^7 ( - , — , . 



Ja 



r 



p C5 

7?' (a2T)') 



+ 7 



Finally, for terms of type (3) we have 



yPg— "g-C6(«/c7)''g-7(c7/w) 



du 



< 



j|7|c7a"Tg-C7/(a"T) 
Ro(p) 



^Re(p)g-C6(«/C7)'' 



C7/(a2T) 



Thus 



lim 

7-i-O 



„ Re(p)/7/ 



A similar calculation yields 



lim 

7-i-O 



dt 



< 



-r 



1 f^eip) C6 

< — „ , , , 7 ' 



After putting in the numbers we obtain 



where 



+h{a,b), 



1/4(0, 6) I < — — 



1 C 



2 aC('*™-i)/2 V 2 'TV(d»-i) 
and lim I^{a,b) = 0. This concludes the proof. □ 

a— ^00 



+ 4 (a, 
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The Casimir pressure is defined to be the force per unit spectral area on each 
parallel plate as the separation b of the two plates is varied, that is, 

The following results are corollaries of the theorem. 
Proposition 7.2. Denote 

As a — > oo, the Casimir energy in a SC waveguide Qa.b is 

\ (-^ 
2^""'" I 2' 



= -4^ 5..(a5CHi E^.^"^'"^^^ +^("^^)- (7-5) 

T/ie zero-temperature Casimir pressure on each parallel plate a{GSC) x {0} or 
a{GSC) X {b} is 



(47r) 2 Geo pez 



log it 



+0 I a^^ ■iu.(Gsc) ] . (7.6) 



The key take-away message is that the Casimir pressure scales with the sep- 
aration of the plates b as 

Pcas(a,fc) - 6-W'(^^)+i)i/2 (log(^)) as oo, (7.7) 

where i?2(-) is 5 log i?-periodic. We do not know the sign of Pcas as it requires 
more information about the function Gq than currently available. 



7.3. Casimir pressure at positive temperature. To calculate the Casimir pressure 
due to thermal effects, we first derive the partition function, then find the free 
energy and take appropriate derivatives. Since most of the arguments resem- 
ble those in the blackbody radiation calculation (Section |6]), we will omit the 
computational details. 

Lemma 7.3. The .simple poles of the spectral zeta Junction C,q^ ^i'^^) ^'"e con- 
tained in the set 

LJU({^|U{%^|)> (7-8) 

where the residues at the two sets of poles read, respectively, 

Vi^^/d^y 2" r{dk,p/2)- ^'-^^ 
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Proof. It is straightforward to check that only the polynomial terms in the ex- 



pansion of if 6 (0 ( 7.1 1 produce singularities. The rest of the proof is as outlined 
in Theorem |3. 1| and (JoroUary |3.2[ □ 



Lemma 7.4. The thermal terms in the grand canonical partition function 

b,0 /""^ photon gas in the Sierpinski carpet waveguide f2a,b are given by 



log^^„„, = ^J2}2^^y^^^Cid,,, + 2)G,.^, 



2 



where (^X]fc=o Spez) means the double sum excluding the term {k,p) — (d, 0). 
Proof. The arguments are identical to those in the proof of Lemma |6.2| □ 
The Casimir pressure is then defined to be 

19 1 9 „ 

Proposition 7.5. As a oo, the Casimir pressure on each parallel plate 
a{GSC) X {0} or a{GSC) x {5} at inverse temperature /3 reads 

Pc.sA'^,b) = r^''^^"^+'^H, (^™log(^^^^ +o(l), (7.11) 
where H3{-) is ^ log R-periodic with 

""^^"^^T^^^l-/ "^U + ^ + logi^jH ^^logi?j 

(7.12) 



Notice that the dominant term in the positive-temperature Casimir pressure 
is independent of the separation of plates b: like in blackbody radiation, the 
scaling is with respect to temperature. 

This concludes our discussions of various results regarding ideal massive or 
massless Bose gas in Sierpinski carpets. We stress once again the power of the 
heat kernel and the zeta function technique in quantum statistical mechanics. 
Not only does it apply to very general spaces, it also give unambiguous answers 
to major thermodynamic questions, such as the emergence of the Casimir energy. 
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8. Towards interacting Bose gas on Sierpinski carpet graphs 

The above discussion of the ideal massive Bose gas is just the prelude to a full 
understanding of quantum many-body physics on fractals. This is because some 
of the most salient features of BEG (in M.'^ or Z'^), such as superfluidity and 
soliton excitations, cannot be explained by the ideal Bose gas. So it behooves us 
to take the next step and consider models of interacting Bose gas on Sierpinski 
carpets. Mathematically speaking, proving BEG in an interacting Bose system 
is a challenging problem (see [24| for a careful exposition). In this section we 
introduce a particular interacting system which has been studied extensively on 
Z"*, then focus on the so-called hardcore Bose gas, which we hope to study in 
more depth on Sierpinski carpet graphs. 



8.1. The Bose-Hubbard model. In what follows we have in mind a connected 
graph G = [V^E) with bounded degree. To describe an ideal Bose gas on 
G, we use the second-quantized Hamiltonian in the grand canonical ensemble 
dr{—A — fJ.1), where A is the combinatorial graph Laplacian on G. To ex- 
press this Hamiltonian, let aj. and be the boson creation and annihilation 
operators at x e V^, which satisfy the canonical commutation relations (GGR): 
axO^y — a^yttx — Sxy and axCiy — ayUx = for all x,y £ V. We also let rix = a^a^,, 
which has the interpretation as the particle number at x. Then one may verify 
that 



dr{-A - /il) 



{xy)eE 



{alay + alax) - X! ('^ ~ deg(a;))n:r- 



A slight variant of this Hamiltonian uses the adjacency operator in place of 
the graph Laplacian, i.e., 



In a realistic experimental setup, multiple atoms tend not to occupy the same 
site X € V due to forbidden overlap of ato mic orbitals. This can be modelled 
by introducing an on-site repulsion term to (8.1). The resultant Hamiltonian is 
that of the Bose-Hubbard model, 



H 



BH 



{a\ay + alax) - ^I'^Ux + U ^ nx{nx - 1), 



■2) 



{xy)(iE 



x£V 



x£V 



with parameters /i £ M and J, J7 > 0. The Hamiltonian (8.2) is widely used by 
condensed matter theorists to model cold atomic systems in optical lattices. In 
the works of BE G on inhoniogeneous graphs 13 14 25|, the authors investigated 
a special case of 



2) with ^ and C/ = 0. 



34 



Joe P. Chen 



8.2. Hardcore Base gas and the XY model. Let us focus on another special case 



of (8.2 1 where the on-site repulsion C/ = oo. In this scenario for any x £V can 
only take the value or 1, that is, particles interact as if they were hard cores. 
The Hilbert space of this hardcore Bose gas is then (C^)®'', which is spanned by 

the tensor product of the basis vectors l^g) '(^l) }' ^^'P'^'^^'^'^^^^S respec- 
tively the occupation or non-occupation of vertex x. Naturally, the creation and 
annihilation operators at x read 

~ 1^0 Oyl ' " 1^1 

Notice that they do not satisfy the CCR anymore: a^aly — a^^a^ = (1 — 2nj.)5xy. 
We can express the Hamiltonian of the hardcore Bose gas as 



Introduce the Pauli matrices 



Then up to an additive constant (|8.3| becomes 



H/jc = - J ^ {SlSl + SlSy) - ^ ^ S"!, (8.4) 

(xy)eE xev 

which is the Hamiltonian of the quantum (ferromagnetic) XY model, with 
/i playing the role of the external magnetic field. Therefore showing that hardcore 
Bose gas condenses on GSC graphs amounts to showing that the corresponding 
quantum XY model has a phase transition. 

In finding a point of attack to this problem, we review results which are known 



to hold on arbitrary graphs. Given any infinite graph G, 11 showed that if the 
simple random walk on G is recurrent, then for any J > 0, the quantum XY 
model at zero field (/i = 0) on G does not exhibit spontaneous magnetization. 
More precisely, if {G„}„ = {(K,i,i?„)}„ exhausts the recurrent graph G, and 
^n,i3.tj. is the Gibbs state on G„ at (/3, /i), then for all /3 > 0, 



1 



hnr hm j^0Jn,p,f_L 




This result generalizes the Mermin- Wagner theorem (27| from translationally 
invariant lattices to arbitrary graphs. Note that fi = corresponds to bosons at 
half-filling (equal number of occupied and unoccupied vertices on average). 

Specializing to the situation on GSC graphs, we deduce that on any GSC 
graph with ds < 2, the quantum XY model at zero field does not undergo a 
phase transition. Under this setting, the hardcore Bose gas at half-filling does 
not Bose condense. Conversely, if condensation occurs on some GSC graph, then 
the graph is transient. 
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The harder probleni is to prove sufficient conditions for which spontaneous 
magnetization, 



lim Hm inf 



1 



n^oo \Vn\ 



> 0, 



or long-range order, 



hm inf 



\Vn\ 




occtirs for large enough /3. Either condition would imply the non-uniqueness of 
KMS states at low temperatures [sl Section 6.2.6]. 

Problem 8.1. Given any GSC graph G, prove or disprove that the transience of 
simple random walk on G is a sufficient condition for the quantum XY model on 
G to exhibit spontaneous magnetization or long-range order at positive temper- 
atures. In other words, determine whether the hardcore Bose gas at half-filling 
condenses on any GSC graph with ds > 2. 

Recall that on Z'' this has been answered in the positive by [12] . However the 
use of Fourier transform was essential in their proof, which made the generaliza- 
tion to arbitrary graphs difficult. 

Though outside the scope of quantum mechanics, we also mention the abelian 
version of the XY model on general graphs, which is better understood. It has the 
same Hamiltonian (8.4|, but with — {S^, S"^, S^) taking values in the unit 2- 



sphere, i.e., X]j=i('^k)^ = 1- is known that there is no spontaneous rotational 
symmetry breaking at any temperature for recurrent graphs [11', 26 . On the fiip 
side, it appears that transience of simple random walk is a necessary, but not 
sufficient, condition for phase transition in the classical XY model. Y. Peres 
conjectured that the sufficient condition is that there exists a. p G (0, 1) such 
that upon performing an independent bond percolation with pr oba bility p, the 
resultant infinite cluster supports transient simple random walk [28[ Conjecture 
1-9]. 

For GSC graphs, which possess strong symmetry and connectivity, we suspect 
that transience alone is a sufficient condition, though we are not aware of any 
previous work on this matter. The key may be to prove a version of spin-wave 
condensation in the spirit of 15 . 



Problem 8.2. Given any GSC graph G, prove or disprove that the transience 
of simple random walk on G is a sufficient condition for the classical XY model 
on G to exhibit spontaneous magnetization or long-range order at positive tem- 
peratures. 
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